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A concept of a family of unique backward-wave photonic devices, such as frequency up and 
down converting nonlinear-optical mirrors, sensors, modulators, filters and amplifiers is proposed. 
Novel materials are considered, which support coexistence of ordinary and backward waves and 
thus enable enhanced nonlinear-optical frequency-conversion processes. Particular properties of 
short-pulse regime are investigated. 



I. INTRODUCTION 

Optical electromagnetic radiation plays important role 
in pollution, chemical and image sensing. Hence, im- 
provement of optical sensors is important for the indi- 
cated branches of sensing. This work is to further develop 
a concept of ail-optically controlled, remotely actuated 
and interrogated sensor that can be employed for envi- 
ronmental probing in remote or hostile locations |l| . Ex- 
traordinary frequency-conversion processes in novel type 
of nonlinear-optical (NLO) materials are proposed and in- 
vestigated that enable great enhancement and frequency 
conversion of signals, which carry important information, 
to the range of maximum sensitivity of detectors. Concur- 
rently, the signals can be redirected towards the detectors. 
Particularly, extraordinary properties of coherent NLO 
energy exchange between ordinary and backward waves 
(BWs) in short pulse regime are investigated. 

The possibilities of enhanced coherent energy exchange 
between electromagnetic waves and amplification of sig- 
nals originate from backwardness, the extraordinary 
property that electromagnetic waves acquire in nega- 
tive index metamaterials (NIMs). Unlike ordinary waves 
propagating in positive-index materials, the energy flow, 
S, and the wave- vector, k, become counter-directed in 
BWs, which determines their unique linear and NLO 
propagation properties. Usually, NIMs are nanostruc- 
tured metal-insulator composites with a special design 
of their building blocks at the nanoscale. Metal com- 
ponent imposes strong absorption of optical radiation in 
NIMs, which presents a major obstacle towards their nu- 
merous prospective exciting applications. Extraordinary 
features of coherent NLO energy conversion processes in 
NIMs that stem from wave-mixing of ordinary and back- 
ward electromagnetic waves (BEMWs) and the possibil- 
ities to apply them for compensating the outlined losses 
have been predicted (for a review, see @, Q an d references 
therein). Most remarkable feature is distributed feed- 
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back behavior of output amplified and generated beams 
which allows for sharp resonance increase of the conver- 
sion efficiency for the appropriate intensity of fundamen- 
tal beam. Essentially different properties of three-wave 
mixing (TWM) and second harmonic generation have 
been shown. Here, we propose two different classes of such 
materials: metamaterials (MM) with specially engineered 
spatial dispersion of their nanoscopic structural elements 
and crystals that support optical phonons with negative 
group velocity. Both do not rely on nanoresonators which 
provide negative optical magnetism and constitute cur- 
rent mainstream in fabricating NIMs. The appearance 
of BEMWs in metaslabs made of standing carbon nan- 
otubes is shown. Uncommon phenomenon of generating 
of a contra-propagating wave at appreciably different fre- 
quency in the direction of reflection is investigated in the 
first class of MM. The possibility to replace plasmonic 
NLO MMs, which are very challenging to fabricate, by 
the ordinary, readily available crystals is another pro- 
posed option. The possibility to mimic unparallel NLO 
frequency-conversion propagation processes attributed to 
NIMs is shown for some of such crystals whereby opti- 
cal phonons with negative group velocity and a proper 
phase-matching geometry are implemented. Here, opti- 
cal phonons are employed instead of BEMWs. The focus 
is on specific properties of amplification and nonlinear- 
optical energy exchange between the contra-propagating 
short pulses. 



II. BASIC IDEA 

A. Huge Enhancement of Nonlinear Optical Energy 

Conversion and Signal Amplification through 

Three-wave Mixing of Ordinary and Backward 

Electromagnetic Waves 

An exotic electromagnetic property of NIMs stems from 
the fact that energy flow and phase velocity of electro- 
magnetic waves become counter-directed inside the NIM 
slab. Such phenomenon of backwardness of electromag- 
netic waves does not exist in naturally occurring materi- 



als. The appearance of BEMW can be explain as follows. 
The direction of the wave-vector k with respect to the 
energy flow (Poynting vector) depends on the signs of 
electrical permittivity e and magnetic permeability \x: 

S = (c/4tt)[E x H] = (c 2 k/4:TTU}e)H 2 s = (c 2 k/47ru;/i)£ 2 . 

If e < and /i < 0, refractive index becomes negative, 
n = —y/JIi, and vectors S and k become contradirected, 
which is in striking contrast with the electrodynamics of 
ordinary, positive index (PI) media (PIM). Hence, mag- 
netic response at optical frequencies, including magnetic 
nonlinear polarization, opens new avenues in electromag- 
netics and for its numerous revolutionary breakthrough 
applications. Such property does not exist in naturally 
occurring materials but becomes achievable in the plas- 
monic metamaterials. Nonlinear-optical propagation pro- 
cesses intrinsic to mixing of ordinary and backward elec- 
tromagnetic waves have been shown to possess unprece- 
dented properties 0-0|. 

Unusual features of NLO coupling of the ordinary and 
backward EM waves that enable compensating losses for 
NI signal through optical parametric amplification (OPA) 
can be summarized as follows. Usually, plasmonic meta- 
materials possess NI properties only within a certain fre- 
quency band and ordinary, PI properties beyond such 
a band. Figure HJa) depicts a slab of thickness L with 

(2) 

quadratic magnetic NLO response x-m ■ The medium with 
quadratic electric response would exhibit similar behav- 
ior. A NI (signal) wave Hi at u)\ enters the slab at its rear 
interface and a strong PI control field, Hz at u>z, enters the 
slab at its front interface. The control field at W3 and weak 
signal at uj\ then generate a difference-frequency idler at 
uj-2 = CJ3 — wi which falls in the PI frequency domain. 
The idler contributes back to the signal through the sim- 
ilar TWM process, u)\ — cj 3 — uj 2 , and thus provides OPA 
of the signal. Such choice of the propagation directions 
ensures all wave vectors ki^.3 to appear co-directed. 

Crucial importance of the outlined geometrical reso- 
nances and striking difference of NLO propagation pro- 
cesses in negative index double domain NIM compared 
with their counterparts in ordinary materials are illus- 
trated in Fig. Q] adopted from our paper [l(|. Besides the 
factor g, the local NLO energy conversion rate for the sig- 
nal is proportional to the amplitude of the idler (and vice 
versa) and depends on the phase mismatch Afc. Hence, 
the fact that the waves decay in opposite directions causes 
a specific, strong dependence of the entire propagation 
process and, consequently, of the transmission properties 
of the slab on the ratio of the signal and the idler decay 
rates. Such extraordinary resonance behavior, which oc- 
curs due to the backwardness of the light waves in NIMs, 
is explicitly seen when compared with similar distribu- 
tion in ordinary, PI materials depicted in Fig. [IJb). Basi- 
cally, such induced transparency resonances are narrow, 
like those depicted in Fig.QTb) and by the plot in Fig.QTc) 
corresponding to a 2 L = 1. This indicates that the sam- 
ple remains opaque anywhere beyond the resonance val- 
ues of the control field and of the sample thickness. Any 
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FIG. 1. (a) Coupling geometry. Si - negative index signal, 
S3 - positive-index control field, S2 - positive-index idler, (b) 
Solid: the dependence of the output signal T\(z = 0) on the 
slab thickness and on intensity of the control field (on factor 
gL). Here, the metamaterial is absorptive at the frequency 
of the signal and assumed amplifying for the idler, Afc = 0. 
Dashed: the same dependence for the ordinary, PI media with 
the same other medium parameters, (c) Transmission T\(z = 
0) of the signal at ol\L — 2.3 and different values of a^L. 
AkL = 0. (d-f) Effect of phase-mismatch on the output signal. 
ot\L — 2.3, oliL = 3, AkL = n. Diminished detrimental effect 
of phase mismatch on the transmission with the increase of 
intensity of the control field and/or the slab thickness is seen. 



sharp frequency dependence of nonlinear susceptibility or 
absorption indices translates into frequency narrow-band 
filtering. The slab becomes transparent within the broad 
range of the slab thickness and the control field intensity if 
the transmission in all of the minimums is about or more 
than 1. Figures [IJc,d) show the feasibility of achieving 
robust transparency and amplification in a NIM slab at 
the signal frequency through a wide range of the control 
field intensities and slab thicknesses by the appropriate 
adjustment of the absorption indices a-i > ql\. Oscillation 
amplitudes grow sharply near the resonances, which in- 
dicates the possibility of cavity-less self- oscillations. The 
distribution of the signal and the idler inside the slab 
would also dramatically change. These simulations prove 
that, unless optimized, the signal maximum inside the 
slab may appear much greater than its output value at 
z = 0. Giant enhancement in the resonances indicates 
that strong absorption of the left-handed, negative-phase 
wave and of the idler can be turned into transparency, am- 
plification and even into cavity-free self-oscillation. Self- 
oscillations would provide for the generation of entangled 
counter-propagating left-handed, fiwi, and right-handed, 
hui2i photons without a cavity. The outlined features 
can be employed for design of ultracompact optical sen- 



sors, selective filters, amplifiers and oscillators generating 
beams of counterpropagating entangled photons. 



B. Three Alternative Coupling Schemes — Three 
Sensing Options 

Figure [2] depict three possible options for the phase 
matched NLO coupling of the ordinary and backward 
waves. Consider the example depicted in panel (c). As- 
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upconverted waves are sent in the opposite directions to- 
wards the remote detectors. Such a signal can be, e.g., 
infrared thermal radiation emitted by the object of inter- 
est, or signal that carries important spectral information 
about the chemical composition of the environment. The 
schemes depicted in Fig. IHa),(b) offer two other options 
with different operational properties for nonlinear-optical 
sensing. The research challenge is that such unprece- 
dented NLO coupling schemes lead to changes in the set 
of coupled nonlinear propagation equations and boundary 
conditions as compared with the standard ones known 
from the literature. This, in turn, results in dramatic 
changes in their solutions and in multiparameter depen- 
dencies of the operational properties of the proposed sen- 
sor. 



III. COHERENT NONLINEAR OPTICAL 

COUPLING OF ORDINARY AND BACKWARD 

ELECTROMAGNETIC WAVES IN 

METAMATERIAL MADE OF CARBON 

NANOTUBES 



FIG. 2. Three different options of the proposed NLO sensors. 
Si,2 and ki,2 are energy fluxes and wavevectors for the ordi- 
nary, positive index, signal and generated idler; S3 and k3 - 
for the negative index control field, (a) Frequency-converting 
metaamplifier. (b) and (c) frequency-converting and ampli- 
fying metamirrors. (b) The NLO sensor amplifies the signal 
Si travailing against the control beam \n(un) < 0], frequency 
up-converts (S2) and sends it along the control beam Si. (c) 
The NLO sensor amplifies the signal Si traveling along the 
control field, up-converts its frequency of and sends it back in 
the direction against the control beam [n(wa) < 0]. 

sume that the wave at uq with the wave-vector ki di- 
rected along the z-axis is a PI (n\ > 0) signal. Usually 
it experiences strong absorption caused by metal inclu- 
sions. The medium is supposed to possess a quadratic 
nonlinearity x and is illuminated by the strong higher 
frequency control field at 0J3, which also falls into the PI 
domain. Due to the three-wave mixing (TWM) interac- 
tion, the control and the signal fields generate a difference- 
frequency idler at UJ2 = W3 — wj., which is a NI wave 
(ri2 < 0). The idler, in cooperation with the control field, 
contributes back into the wave at uq through the same 
type of TWM interaction and thus enables optical para- 
metric amplification (OPA) at uq and enhanced reflected 
beam at 102 by converting the energy of the control fields 
into the signal and the idler. In order to ensure effective 
energy conversion, the traveling wave of nonlinear polar- 
ization of the medium and the coupled electromagnetic 
wave at the same frequency must be phase matched, i.e., 
must meet the requirement of Ak = k3 — k2 ki =0. 
Hence, all phase velocities (wave vectors) must be co- 
directed. Since nioJi) < 0, the idler is a BW, i.e., its en- 
ergy flow S2 = (c/47r)[E2 x H2] appears directed against 
the z-axis. This allows to conveniently remotely inter- 
rogate the NLO chip and to actuate amplification and 
frequency up-conversion of the signal. Transmitted and 



Exciting unparallel avenues for nonlinear electromag- 
netics can be open by the metamaterials where formation 
of B Ws becomes possible due to specific spatial dispersion 
of their structural elements. Basic idea is as follows. As 
outlined above, according to currently commonly adopted 
concept, negative refractive index and associated back- 
wardness of optical waves require negative permeability 
and therefore magnetism at optical frequencies. How- 
ever, a different approach is possible [ll|, [12] ■ In a loss- 
free isotropic medium, energy flux S is directed along the 
group velocity v g : S = v g U, v g — gradkw(k). Here, U 
is energy density attributed to EMW. It is seen that the 
group velocity may become directed against the wavevec- 
tor depending on sign of dispersion duj/dk. Basically, 
negative dispersion can appear in fully dielectric materi- 
als with particular dispersion of its structural elements. 
This opens an entirely novel research and application av- 
enue. In the next subsection, we show a possibility to 
engineer such a dispersive medium which supports coex- 
istence of both ordinary and BEMWs which can be phase 
matched. An example of enhanced coherent energy ex- 
change between ordinary fundamental EMW and its sec- 
ond harmonic is considered in a lossless NLO slab for the 
sake of simplicity. 



A. Carbon "Nanoforest" and Phase Matching of 
Ordinary and Backward EM Waves. 

Appearance of BEM modes in nanoarrays and layered 
structures has been predicted recently in Ref. |13l - tl5| . 
Obviously, many other options should have existed. Be- 
low, we propose such an option that seems promising in 
the context of nonlinear propagation and coherent energy 
conversion processes. Namely, the possibility of conver- 
sion of ordinary EMW to the contrapropagating BEMW 



at its doubled frequency. Hence, such metaslab can be 
viewed as a frequency-doubling NLO mctamirror. 
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FIG. 3. (a) Geometry of free-standing CNTs. (b) Dispersion 
- the frequency vs slow- wave factor for the slab of CNTs with 
open ends, (c) Group delay factor vs the phase velocity slow- 
wave factor for the same modes as in panel (a). Black (flat) 
curve corresponds to the high-frequency mode, blue curve to 
the low- frequency mode. The tip of blue curve is cut. Its 
maximum corresponds to the stop-light regime. 



y/tzz [k 2 — (rmr/2h) 2 ], where m is a positive integer, h is 
the height of the waveguide (CNT) and k is the wavenum- 
ber in free space. If e zz < 0, BW propagation is al- 
lowed when k < rmr/2h and forbidden for k > mir/2h. 
The relation between the wavevector component k x and 
wavenumber k is: k 2 = \(k 2 — k p )(k 2 — k 2 )]/[k 2 ], where 
k z = rmr/(2h), m is the integer determining a number of 
field variations along CNT, k p is plasma wavevector. One 
can show that dfc^/dfc 2 < 0, if k z /k > 1 and k p /k > 1. 

Numerical analysis of Eq. ([2]) is depicted in Fig. [3jb) 
for the case of CNT radius r = 0.82 nm, the lattice 
period d = 15 nm and EM modes with m = 1 and 
m = 3. Slow-wave factor is proportional to the wave 
vector: c/v p h = k(u)/ko, where ko is the wave vector in 
the vacuum. The appearance of positive dispersion for 
small slow-wave factors is caused by interaction of BW 
in the CNT slab with the plane wave in air. Indeed, 
coexistence of the positive (ascending dependence) and 
negative (descending dependence) dispersion for differ- 
ent frequencies proves that such a metamaterial supports 
both ordinary and backward EMW. It also proves that 
resonant plasmonic structures, like split-ring resonators, 
exhibiting negative e and fj, are not the necessary re- 
quirement for the realization of BW regime in mid-IR 
range. The possibility of considerable increased band- 
width of BEMW compared to most plasmonic MM made 
of nanoscopic resonators is seen that gives the ground 
to consider CNT arrays as a promising perfect backward- 
wave metamaterial. The slow-wave factor for both modes 
is shown in Fig. OUc). The magnitude of n gr goes to infin- 
ity at n p h rj 1.85, which indicates the stop-light regime 
for the low-frequency mode. Particularly, Fig.[3Jb) shows 
the possibility of phase matching of ordinary fundamental 
and backward second harmonic EMWs. 



Figure [3ja) depicts a periodic array of carbon nan- 
otubes (CNT) vertically standing on the surface of a per- 
fect electric conductor (PEC) with the CNT ends open 
to air, which can be seen as per fect magnetic conduc- 
tor (PMC). As shown in Ref. [15|], EM waves travelling 
through such CNT "nanoforest," along x or y directions, 
posses a hyperbolic dispersion and relatively low losses in 
the THz and mid-IR ranges. One of the most important 
consequences from the hyperbolic-type dispersion law is 
the possibility for propagation of both forward and back- 
ward EM waves. Consider EMW propagating along the 
x-axis. We also introduce group n gr — c/v g and phase 
n p h = c/v p h slow-wave factors. The latter one is refrac- 
tive index. For the given case of surface waves propa- 
gating in the slab of CNTs with open ends, whose fields 
attenuate in air, the dispersion is given by the equation: 



tsai(k z h) = \/k^ — k 2 /k z . 



(2) 



Such a dependence can be understood from considering 
a planar waveguide formed by perfect PEC and PMC 
planes and tampered with a CNT array. The array 
axis is orthogonal to the walls of the waveguide. Then, 
the propagation constant along the waveguide is k± = 



B. Coherent Energy Exchange Between Short 
Contrapropagating Pulses in the Carbon Nanoforest 

Here, we demonstrate unusual dependence of pulse 
shape and overall efficiency of SHG on the ratio of input 
fundamental pulse length and metamaterial slab thick- 
ness. Consider a double domain positive/negative index 
slab of thickness L that supports ordinary EMW at fun- 
damental frequency (FH) and BEMW at SH frequency. 
To ensure phase matching, wavevectors of the FH and 
SH must be co-directed and, hence, their energy fluxes 
- counter-directed. The slab operates as a frequency up- 
converting nonlinear- optical mirror with controllable re- 
flectivity. The equations for amplitudes of FH, <zi, and 
SH, a2, can be written as 



-i2ga*a2 exp (iAkz) - 
= igaf exp (—iAkz) 



— ( 
2 

Oil 



-a 2 . 



(3) 
(4) 



1 da\ da\ 
vi dt dz 

1 da 2 da 2 
V2 dt dz 

Here, |ai,2| 2 are slowly varying amplitudes proportional to 
the instant photon numbers in the energy fluxes, ai.2 are 
absorption indices, Ak = k 2 — 2ki is phase mismatch, and 



Vi are the group velocities for the corresponding pulses. 
Note opposite signs in the equations and the requirement 
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FIG. 4. Input, Ti(z = 0), and output, Ti(z — L), pulses 
of fundamental and negative- index SH radiation, 772(2 = 0). 
(a) and (b) - short pulse: d=4. (c) and (d): longer pulse: 
d=l. (b) and (d) - input power is 25 times larger than in (a) 
and (c). (a) Input pulse area (energy) S10 = 0.9750; output 
pulse areas (energy) Sil — 0.5031, S20 = 0.2392. (b) Input 
pulse energy S10 = 0.9750; output pulse energy Sil — 0.0396, 
S20 = 0.4742. (c) Input pulse energy S10 - 0.9900; output 
pulse energy Sil = 0.2516, S20 = 0.3692. (d) Input pulse 
energy S10 = 0.9900; output pulse energy Sil = 0.0161, $20 = 
0.4870. 

of the boundary conditions to be set at the opposite edges 
of the slab for FH and SH. These lead to cardinal changes 
in the solutions to the equations as compared to those in 
the ordinary nonlinear optical material. We have chosen 
the input pulse shape as being close to a rectangular form 



F(t) = 0.5 [ tanh -^-±-1 
\ St 



tanh 



to 



St 



(5) 



where St is the duration of the front and tail, and ro is 
the shift of the front relative to t = 0. All quantities 
are reduced by the pulse duration At. The magnitudes 
St = 0.01 and To = 0.5 have been selected for numerical 
simulations. 

Unusual properties of BWSHG in NIMs in the pulsed 
regime stem from the fact that it occurs only inside the 
traveling pulse of FH. SHG begins on its leading edge, 
grows towards the back edge, and then exits the pulse 
with no further changes. Since the FH pulse propagates 
across the slab, the duration of the SH pulse is longer 
than the fundamental one. Depletion of the FH radiation 
along the pulse length and the conversion efficiency de- 



pend on its initial maximum intensity and phase match- 
ing. Ultimately, the overall properties of SHG, such as 
the pulse length, and the photon conversion efficiency, 
appear dependent on the ratio of the FH pulse and slab 
lengths. Such an extraordinary behavior is illustrated in 
Figs. Ufa)-(d). Here, d is the slub thickness reduced by 
the input pulse length, d = L/v\At, g is proportional 
to product of nonlinear susceptibility x^ 2 ) and the in- 
put amplitude of FH. A rectangular shape of the input 
FH pulse T\ = |ai(z)| 2 /|aio| 2 is depicted at z = when 
its leading front enters the medium. The results of nu- 
merical simulations for the output FH pulse, when its 
tail passes the slab's edge at z = L, as well as for the 
shape and conversion efficiency of the output SH pulse, 
V2 = |a2(z)| 2 /|aio| 2 , when its tail passes the slab's edge 
at z = 0, are shown. For clarity, here, the medium is as- 
sumed loss-free, group velocities of the fundamental and 
SH pulses assumed equal, Afc = 0. Panels (a) and (b) cor- 
respond to the fundamental pulse four time shorter than 
the slab thickness. Increase of the conversion efficiency 
with increase of the intensity of the input pulse is seen. 
It is followed by shortening of the SH pulse. The outlined 
properties satisfy to the conservation law: the number of 
annihilated pair of photons of FH radiation (S10 — Sxl)/2 
is equal to the number of output SH photons S^o- Pan- 
els (c) and (d) display corresponding changes for a longer 
input pulse with the length equal to the slab thickness. 
Here, larger conversion efficiency can be achieved at a 
lower input intensity compared with the preceding case 
because of the longer conversion length. The changes in 
the SH pulse length and conversion efficiency with in- 
crease of input intensity appear less significant. 



IV. FULLY DIELECTRIC BACKWARD- WAVE 
NLO MATERIAL: ENHANCING COHERENT 

ENERGY TRANSFER BETWEEN 

ELECTROMAGNETIC WAVES IN ORDINARY 

CRYSTALS BY COUPLING WITH OPTICAL 

PHONONS WITH NEGATIVE PHASE VELOCITY 

As above described, NLO with backward electromag- 
netic waves enables a great enhancement of energy- 
conversion rate at the otherwise equal nonlinearities and 
intensities of input waves. Herein, we propose fundamen- 
tally different scheme of TWM of ordinary and backward 
waves (BW). It builds on the stimulated Raman scat- 
tering (SRS) where two ordinary EM waves excite back- 
ward elastic vibrational wave in a crystal, which results 
in TWM. The possibility of such BWs was predicted by 
L. I. Mandelstam in 1945 [lf|, who also had pointed out 
that negative refraction is a general property of the BWs. 
The idea underlying the proposed concept and its ba- 
sic justification is described below. The goal is to show 
the possibility to replace the NI plasmonic composites, 
which are challenging to fabricate, with readily available 
ordinary crystals, some of which have been already exten- 
sively studied, and thus to mimic the unparallel properties 
of coherent NLO energy exchange between the ordinary 



and BW. The basic idea is as follows. In 17], stimulated 
Raman scattering (SRS) on optical phonon was investi- 
gated in continuous wave (CW) regime. The possibility of 
distributcd-fcedback type resonance enhancement of am- 
plification has been shown that stem from negative disper- 
sion of elastic waves. The effect appeared similar to that 
in three- wave mixing (TWM) of ordinary and BEMWs in 
NIMs, provided that special requirements are met. Yet, 
the required intensity of the fundamental field was found 
to be close to the optical breakdown threshold due to the 
high phonon damping rate. Here, we show that the in- 
dicated fundamental formidable obstacle can be removed 
by making use of short pulses, which opens the possibility 
to mimic TWM processes attributed to NIMs in readily 
available crystals. Besides that, we show that such mixing 
exhibits unparallel properties that allow for tailoring the 
shapes of the generated and transmitted pulses and for 
huge enhancement of the frequency conversion efficiency. 




FIG. 5. Negative dispersion of optical phonons and two phase 
matching options for short- and long-wave vibrations: (a) - co- 
propagating, (b) - contra-propagating fundamental (control) 
and Stokes (signal) waves. Insets: relative directions of the 
energy flows and the wave-vectors. 

Typical dispersion curve cj(fc) for optical phonons, 
which exist in crystals containing more than one atom 
per unit cell, is depicted in Fig. [5] The dispersion 
is negative in the range from zero to the boundary of 
the first Brillouin's zone. Hence, the group velocity of 
such phonons v„ is antiparallel with respect to its wave- 
vector h v . The dispersion w v (k v ) can be approximated as 
oj v = y/uiQ — 0k% . Then, in the vicinity of k v = 0, veloc- 
ity v^ r is given by Ujj r = —(3k v /uj v — —(3/vP h , where v% h 
is the projection of the phase velocity of the vibrational 
wave on the z-axis and /3 is the dispersion parameter for 
the given crystal, eql Optical elastic vibrations can be ex- 
cited by the light waves through the Raman scattering. 
The latter gives the ground to consider such a crystal as 
the analog of the medium with negative refractive index 
at the phonon frequency and to employ the processes of 
parametric interaction of three waves, two of which are or- 
dinary EM waves and the third one is the backward wave 
of elastic vibrations. The coupled waves are described by 
the equations 

E li8 = (l/2)£ liS (z, t)e ikl -''- iUl " t + c.c, (6) 

Q v = (l/2)Q(z, tje**-*"*"-* + c.c. (7) 

Here, £^ s , Q, loi sv and ki <StV are the amplitudes, frequen- 
cies and wave-vectors of the fundamental, Stokes and vi- 
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FIG. 6. (a) Transmission of the Stokes wave T^ (z = L) vs 
intensity of the fundamental control field in the vicinity of first 
"geometrical" resonance [co-propagating Ei and E s geometry, 
Fig. Ufa)] . g = \/gZg s . Such extraordinary resonance appears 
because of backwardness of the coupled vibration wave and 
opposite direction of propagation for the Stokes and phonon 
waves, (b) Comparison of the output intensities of the Stokes 
wave vs intensity of the control field for co- (the blue, solid 
line) and contra-propagating (the red, dashed line) fundamen- 
tal (control) and signal (Stokes) waves [coupling geometries 
depicted in Fig. [SJa) and (b) respectively] . 



brational waves; Q v (z,t) — yfpx(z,t); x is displacement 
of the vibrating particles, p is the medium density and the 
requirements uji = u> s + uj v (k v ), ki = k s + k v are sup- 
posed met. Partial differential equations for the slowly 
varying amplitudes in the approximation of the first or- 
der of Q in the polarization expansion are [18] : 
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Here, vi, s and —v v are projections of the group velocities 
of the fundamental, Stokes and vibration waves on the z- 
axis, N is the number density of the vibrating molecules, 
a is the molecule polarizability, t v is phonon lifetime. 

Equations (|8|)-(fT0|) are similar to those describing 
TWM of contra-propagating waves in NIMs for the phase- 
matching scheme of co-propagating fundamental and 
Stokes waves, where vibration wave counter-directed to 
the Stokes wave, which is depicted in Fig. 03a). In the 
case of continuous waves and neglected depletion of the 
fundamental wave, resonance enhancement becomes pos- 
sible, similar to those, depicted in Fig. [1] . On the con- 
trary, only standard exponential behavior is possible in 
the case of Fig. [5fb). Corresponding transmission fac- 
tors in the vicinity of first "geometrical" resonance are 
shown in Fig. [5] In the resonance, Tj* —¥ oo, which is 
due to the approximation of constant control field. Con- 
version of the control field to the Stokes one and excited 
molecule vibrations would lead to saturation of the con- 
trol field which limits the maximum achievable amplifica- 
tion. Strong amplification in the maximums indicates the 
possibility of self-oscillations and thus creation of mir- 
rorless optical parametrical oscillator with unparalleled 
properties. Figures H{a,b) indicate the possibility to fit 



in the effective conversion length within the crystal of a 
given thickness and to significantly concentrate the gen- 
erated Stokes field nearby its output facet. Such atypical 
extraordinary behavior in readily available crystals may 
find exciting applications. However, the estimates have 
shown that intensity of the fundamental field which is to 
attain such extraordinary amplification appears close to 
the optical breakdown threshold [17j . It is because of fast 
phonon damping and corresponding high rate of energy 
conversion of the fundamental beam in heat. 

Below, we show that such seemingly unavoidable ob- 
stacles can be overcome in short pulse regime. For the 
sake of simplicity, we consider model of a rectangular 
pulse of input fundamental radiation with the pulse du- 
ration t p <C t v . A seeding Stokes wave is assumed a 
weak continuous wave. In the moving coordinate frame 
associated with this pulse and within its range, complex 
amplitudes of two other interacting fields become time in- 
dependent. Then the analytical solution to the equations 
can be found for amplification corresponding to negligible 
depletion of the pump due to the conversion. Numerical 
solution is found for the opposite case. Inside the crys- 
tal, the boundary conditions must be fulfilled not at the 
boundaries of the medium but at the boundaries of the 
fundamental pulse. The latter is correct for the period 
of time after the instant when the generated Stokes and 
vibration pulses reach the boundaries of the fundamental 
pulse. Such approximation becomes true after travail- 
ing a distance I > i max , where L max = max{L s ,L v }, 
L s = l p vi/\v B - vi\ and L s > v = l p vi/(vi + v v ), l p = T p v t is 
length of the fundamental pulse. Hereinafter, the waves 
are referred to as co-propagating if the Poynting vector 
of the Stokes wave is co-directed with that of fundamen- 
tal wave [Fig. [SJa)], and as counter-propagating in the 
opposite case [Fig. Ejjb)]. Note, that v s is negative for 
contra-propagating Stokes wave. 

In the approximation of constant pump amplitude, 
in the coordinate frame locked to the pump pulse, and 
within the pulse, equations for the generated Stokes and 
backward vibration waves take the form: 

dQ/dt; = —ig v £* + QK v /l v , dE s /dS, = ig s Q* ■ (11) 

Here, £ = z—vft, l v = t v v v is the phonon mean free path, 
g v = K v N(da/dQ)£i/(Auj v v v ), K v = v v /(vi +v v ), g s = 
K s N(da/dQ)£inuj 2 s /(k s c 2 ), K s = v s /(v s - v t ); v s ,k s > 
for co-propagating and v s , k s < for counter-propagating 
beams. Since the Stokes frequency is less than that of the 
fundamental one, v s > vi and vi 3> v v . 

Equations (TTT|) are similar to those describing CW 
TWM, [17| except the boundary conditions. They cor- 
rectly describe possible huge amplification of the Stokes 
signal until relatively small part of the strong input laser 
beam is converted. For co-directed laser and Stokes 
waves, the boundary conditions are: 



£ s (£ = 0)=£ s °, Q(£ = / p )=0. 



(12) 



In the opposite case, they are written as 

£.(£ = *„)=#, Q(Z = l P ) = 0. (13) 



The analysis of solution to Eqs. (jlip shows that, in the 
given approximation of neglected depletion of the funda- 
mental wave, amplification of co- directed signal tends to 
infinity, when the pulse energy approaches the resonance 
value corresponding to gl p = tt/2, where g = y/g^gs- This 
indicates the possibility of huge enhancement of the con- 
version efficiency. Respective intensity of the fundamental 
field P- is given by the equation 



K v cn s X s0 ui v 



mm K s 16tt3^t2 



tp - 



N 



da 
dQ 



(14) 



where, n s is refractive index at uj v and A s o is the wave- 
length in vacuum. From comparison with the correspond- 
ing value I m in for the CW regime [17]], one concludes 
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For the cr ysta l parameters, which are characteristic 
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for calcite [19( and diamond [2flH22|, Eq. (JEJ) yields 
^min/^min pa 10~ u and, hence, suggests a decreases of 
J min down t0 J min ~ l ° 7 W/cm 2 . The latter is achievable 
with commercial femtosecond lasers and falls below the 
optical breakdown threshold for most transparent crys- 
tals. 

Equation (TIB]) displays two factors that determine sub- 
stantial decrease of I^ in in pulsed regime compared to 
that in CW one. First factor is the ratio of group velocity 
of the elastic wave to that of the fundamental one, v v /vi, 
which is on the order of ~ 1CP 8 . This factor is attributed 
to the fact that phonons generated on the front edge of 
the laser pulse propagate in the opposite direction and, 
hence, exit very fast, practically with the optical group 
velocity, from the fundamental pulse zone before it is dis- 
sipated. Hence, effective phonon mean free pass becomes 
commensurable with the fundamental pulse length. This 
mitigates the detrimental effect of phonon damping. With 
increase of v v , phonon mean free path grows, which de- 
creases both I^ in and / m ; n in a way that the advantage 
of pulse regime over CW regime diminishes. The second 
factor in Eq. (fTS")) determines further decrease of i^ in due 
to small optical dispersion in the transparency region of 
the crystals. The fact that the Stokes pulse surpasses the 
fundamental one slowly increases significantly the effec- 
tive NLO coupling length. 

To investigate the regimes of significant energy conver- 
sion, a set of partial differential Eqs. (I5))- (|10[) was solved 
numerically in three steps: TWM in the vicinity of the 
entrance, inside and in the vicinity of the exit from the 
Raman slab. Simulations for the first and third inter- 
vals were made in the laboratory reference frame with 
the boundary conditions applied to the corresponding 
edges of the slab. The propagation process inside the slab 
was simulated in the moving frame of reference with the 
boundary conditions applied to the pulse edges. Such an 
approach allowed for significant reduction of the computa- 
tion time because, for each given instant, the integration 
was required only through the space interval inside the 
fundamental pulse and not through the entire medium. 
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FIG. 7. a. Conversion efficiency vs. intensity of the input 
pump for co-propagating (solid line)ss and contra-propagating 
(dash line) geometries, b. Output Stokes (solid line) and 
fundamental (dash line) pulses for co-propagating coupling. 



Shape of the fundamental pulse was chosen nearly rect- 
angular and symmetric with respect to its center 

Si = -£;°{tanh[(t +t p - t)/t f ] - tanh[(t - *)/*/]}• 

The slope tf — 0.1, the pulse duration at half-maximum 
t p = 1 and its delay to = 0.6 were scaled to the fundamen- 
tal pulse width r p . The amplitude of input CW Stokes 
signal was chosen £® = 10~ 5 £ ; °. Numerical investiga- 
tions were done for the model with parameters typical 
for calcite [19] {lo v — 1086 cm -1 ) and diamond [20l423 | 
(lu v = 1332 cm -1 ): carrying wavelength A; = 800 nm, 
pulse duration t p — 60 fs, u) v = 1332 cm -1 , t v — 7ps, vi = 
1.228 • 10 10 cm/s, v s = 1.234 • 10 10 cm/s, v v = 100 cm/s 
for co-propagating and v v = 2000 cm/s for counter- 
propagating waves, Nda/dQ = 3.78 • 10 7 (g/cm) 1 / 2 , the 
crystal length L=lcm. 

Figure [TJa) displays the output quantum conversion ef- 
ficiency rj q = (u>i/w s ) ■ L I s {z, t)dt/ J t Ii(z — 0, t)dt vs. in- 
put pulse intensity, both for co-propagating (z = L) and 
counter-propagating (z = 0) geometries. Many orders in- 
crease of the conversion efficiency due to BW effect in the 
case of co-propagating waves is explicitly seen. Saturation 
at Ii/I^ in > 7 • 10 4 is due to depletion of fundamental ra- 
diation caused by conversion to Stokes radiation. Figure 
W[b) displays a shortened Stokes output pulse that has 
surpassed the inhomogeneously depleted and broadened 
fundamental pulse (zoomed in the inset). 

The simulations also show that shapes of the output 
Stokes and fundamental pulses differ and vary signifi- 
cantly depending on the intensity of the input funda- 
mental wave. Figure [SJa) depicts amplified output pulse 

T s = \£ s (L,t)/£®\ of Stokes radiation for relatively small 
depletion due to conversion of the fundamental beam. 
Here, shape of the fundamental pulse is unchanged. The 
shape of the amplified Stokes pulse is different and de- 
termined by the fact that v s > vi and Stokes pulse has 
surpassed the fundamental one. In contrast, the quantum 
conversion becomes significant in Figs. IHb)-(d). Corre- 
sponding depletion of the output fundamental pulse and 
changes in its shape are explicitly seen. Note that in the 
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FIG. 8. Changes in the shapes of generated Stokes (solid 
line) and transmitted fundamental (dash line) output co- 
propagating pulses with the increase of energy of the input 
fundamental pulse. r\ q is corresponding conversion efficiency. 



case depicted in Fig. EJd), the output Stokes pulse signif- 
icantly overtakes the pump pulse. In the latter case, ma- 
jor conversion occurred inside, far from the crystal edges, 
and then both pulses propagated without interaction. It 
is seen that output Stokes pulse narrows with the increase 
of energy of the fundamental pulse. Here, crystal length 
of 1cm fits L/lp = 1357 input pulse lengths. The thresh- 
old intensity i^ n = 6 • 10 6 W/cm 2 corresponds to 60 fs 
pulse of about 5 /zJ energy focused to the spot of diameter 
D = 100 /im. Intensity of the seeding Stokes signal was 
chosen / 6 °/7,° = lO -10 . 

Note, that the described NLO propagation process is 
in a striking contrast with their positive group veloc- 
ity counterparts [2j, |24[, including acoustic waves with 
energy fluxes directed against that of EM waves [25j . 
The proposed here concept is different from the one ear- 
lier proposed in [26j and does not require periodic pol- 
ing of quadratic nonlinear susceptibility of crystals at 
the nanoscale as described in [27j (and in the references 
therein) . 



V. CONCLUSION 

The possibility of creation of remotely all-optical con- 
trolled frequency up-converting and amplifying nonlinear 
optical sensor (image and data processing chips) is investi- 
gated and proved with numerical simulations. Novel pho- 
tonic materials are proposed which support coexistence of 



ordinary and backward waves that can be phase matched. 
Such nonlinear-optical materials enable greatly enhanced 
nonlinear-optical frequency-conversion processes that can 
be employed for creation of optical sensors with unparal- 
lel operational properties. Two different classes of such 
materials are proposed: metamaterials with specially en- 
gineered spatial dispersion and crystals that support opti- 
cal phonons with negative phase velocity. Unlike current 
mainstream in fabricating negative index metamaterials, 
both options do not rely on nanoresonators that provide 
negative optical magnetism and thus constitute current 
mainstream in fabricating negative index metamaterials. 
Extraordinary properties of coherent frequency conver- 
sion processes in the proposed materials are numerically 
simulated both in continuous wave and in short pulse 
regimes. 

As an example of dispersion engineering, a metama- 
terial made of standing carbon nanotubes is described 
which supports coexistence of ordinary fundamental and 
backward second harmonic electromagnetic waves that 
can be phase matched. Extraordinary properties of such 
frequency-doubling metamirror operating in short-pulse 
regime appear in striking contrasts with second harmonic 
generation in ordinary materials. 

Fabrication of specially shaped nanostructures which 
enable negative optical magnetism is challenging task 
that relies on sophisticated methods of nanotechnology. 
Engineering a strong fast quadratic NLO response by 
such mesoatoms also presents a challenging goal not 
yet achieved. This work proposes a different paradigm, 
which is to mimic similar extraordinary backward-wave 
frequency conversion processes through making use of 
readily available Raman active crystals. Basic underpin- 
ning idea is to replace one of the coupled backward opti- 
cal electromagnetic mode, which existence is commonly 
attributed to negative-index metamaterials, by optical 
phonons - the elastic wave with negative group velocity. 



Operation in short pulse regime is proposed to remove a 
severe detrimental factor imposed by fast phonon damp- 
ing. Significant decrease of the required minimum inten- 
sity of fundamental radiation, as compared with that in 
the continuous- wave regime, is predicted down to that 
provided by commercial lasers. Unparallel properties of 
the proposed short-pulse process are numerically simu- 
lated and the possibility of huge enhancement of quan- 
tum conversion efficiency as well as of tailoring duration 
and shapes of the generated and transmitted fundamental 
pulses are predicted. 

Further elaboration of the proposed concept allows to 
utilize the revealed extraordinary features for creation 
of a family of unique photonic devices with advanced 
functional properties such as nonlinear-optical mirrors 
with frequency-dependent controllable reflectivity, unidi- 
rectional optical amplifiers, frequency narrow band filters, 
switches and cavity-free optical parametric oscillators, 
which can be exploited for optical sensing. There exist 
many dispersive materials and readily available Raman- 
active transparent crystals that support electromagnetic 
and elastic waves with negative group velocity which can 
be utilized for the given purpose. Some of them may 
offer the opportunities for further optimization of unpar- 
allel operational characteristics of the proposed photonic 
devices. 
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